Spin monopoles with Bose— Einstein condenstates 
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We study the feasibility of preparing a Bose-Einstein con- 
densed sample of atoms in a 2D spin monopole. In this state, 
the atomic internal spins lie in the x-y plane, and point in 
the radial direction. 

PACS: 03.75.Fi,05.30.Jp,32.80.Pj, 



After the successful generation of Bose-Einstein con- 
densation of alkali atoms, the creation and manipulation 
of certain macroscopic quantum states remains as one of 
the fundamental goals in the field of Atomic Physics |Q . 
During the last year, several ways to create vortices and 
solitons have been proposed pLpfl . Under certain circum- 
stances these states are stable |4p|| , which has motivated 
several experimental groups to try corresponding experi- 
ments. In this letter we study a new kind of macroscopic 
quantum state for these atomic samples, what we call a 
2D spin monopole. It is a state in which the atomic spin 
(i.e., the magnetization vector) points in the radial di- 
rection in the x-y plane [Fig. 1(a)]. We show that this 
state is stable under realistic conditions, and analyze a 
method to generate it which only requires current experi- 
mental technology. We will first study a one dimensional 
situation, in which the condensate is confined in a ring, 
and then we will generalize it to the 3D case. 




x,y z 
FIG. 1. (a) Spin monopole in 2D. For fj, ^ the spins rotate 
as a function of time; (b) Ring trap; (c) Potential V(r, z) 

We consider a Bose-Einstein condensed sample of N 
atoms. The atoms have two internal levels, | f) and | |), 
which are relevant for the dynamics (TjJ^]. We assume 
that half of the atoms is in each internal level. The atomic 
motion is confined by an external trap with a ring shape 
[see Fig. 1(b)] |9) . In the limit where the motion along the 



radial (r) and axial (z) direction is frozen, the dynamics 
of the motional state only depends on the polar angle (ff) . 
We can write the wavefunction of the condensate as 



\*(6,t))=M0,t)\ T) + &(0,r)| I), 



(1) 



where 4>i,2 are the motional wavefunctions corresponding 
to the internal states | j, |), respectively and fulfill the 
coupled Gross-Pitaevskii Equations: 



• d ^ 


" d 2 




de 2 




- d 2 




~d9 2 



+ Si+ uii |0i 

+ &2 + 1412 1 </>l 



\ 2 + u 2 M 2 



U22I02I 2 



with normalization 
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Here, #1,2 denotes the energy of the two internal states 
and Uij — Uji (i,j = 1,2) describe their mutual interac- 
tions. All the quantities in Eqs. (]|) have been rescaled 
so that they are dimensionless. In particular, the are 
proportional to the number of atoms TV, the correspond- 
ing scattering length, and the square of the ring radius 
[see Eq. ( pa) ) below]. 
The state (0) with 
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is a stationary solution of Eqs. 



with fij = j — 1 + 



uij + uij + Sj (j = 1,2). Furthermore, defining the Pauli 
operator a — (o~ x , cr y ,a z ) as usual, we have that the state 
([j]) with (|j) is an eigenstate of a ■ n with n = [cos(0 — 
/it), sin(# — fj/r), 0] and /i = fi2 — Hi- This means that for 
fj, = the spins are always pointing outwards in the x-y 
plane, whereas for /i 7^ they oscillate as shown in Fig. 
1(a). Thus, the stationary state (f|) can be considered as 
an oscillating spin monopole in two dimensions. 

In order to analyze the stability of the spin monopole 
we have carried out a perturbative linear analysis. We 
consider a small perturbation around the monopole so- 
lution, such that 411.2(0, t) = (j)^(9,T) + ea\^{0,T). We 
then expand Eq. (^|) up to first order in e obtaining a lin- 
ear set of coupled equations for ai,2 and a* 2 - Expand- 



ing a ia (6,T) = J2n=- 
thesc equations, we obtain 
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is a column vector. Here, TL n = if„ + if 11 
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are 4x4 matrices and <S> denotes tensor product. The 
stability analysis can be carried out by diagonalizing the 
matrices ETL n . Complex eigenvalues ±iA correspond to 
exponentially growing solutions [a(r) = a(0)e AT ], i.e. 
the stationary solution is dynamically unstable. Real 
eigenvalues A lead to small oscillatory solutions [a(r) = 
a(0)e ±lXT ), which may correspond to higher or lower val- 
ues of the Gross-Pitaevskii energy functional with re- 
spect to the stationary solution. In the first case, the 
stationary solution is stable, whereas in the latter one the 
system is energetically unstable; this means that if energy 
can be taken out from the system at a given rate, then 
the stationary solution is unstable on that time scale. 
Another interesting case occurs when one of the matrices 
is not simple (i.e., non diagonalizable), since in that case 
one has solutions which only grow polynomically with 
time [for example, for a 2 x 2 Jordan matrix correspond- 
ing to a degenerate eigenvalue A, a(r) — a(0)(l + At)]. 
On the other hand, if one is only interested in checking 
whether the solution is stable or not, (without distin- 
guishing among the different types of instability) one can 
simply study the positivity of TL n . If all matrices TL n > 
and the projector operator Pq on the kernel of Ti. n com- 
mutes with E for all n, then the stationary solution is 
stable. One can proof this statement as follows: first, 
the positivity of Tl n ensures the positivity of the Gross- 
Pitaevskii energy functional and therefore there cannot 

be energetically instabilities; secondly, all eigenvalues of 

1/2 1/2 

E1-L n are also eigenvalues of Ti n EH. n , which is her- 
mitian (has real eigenvalues) and therefore there cannot 
be dynamical instabilities; finally, if [Pq,E] = one can 
easily prove that if a vector belongs to the kernel (range) 
of Tin, then it also belongs to the kernel (range) of EH n , 
and therefore this last matrix is simple. 

We start out by analyzing the positivity of Ti. n . We 
consider two cases: (a) U11M22 < u 2 2 . The matrix H lnt 
has a negative eigenvalue; for n = 0we have Kq = and 
therefore Hq has a negative eigenvalue and the solution 
is not stable, (b) U11U22 > u i2- Now H lnt > 0. Given 
that K n > for \n\ > 2 we can restrict ourselves to the 
case \n\ — 1. Since the first three diagonal minors of TL n 
are positive, we just have to impose that the determinant 
be positive. We obtain that the monopole ([|) is stable if 
W11M22 — u\ 2 +u 22 /2 — 3iin/2 — 3/4 > 0. Now, we investi- 
gate which kind of instabilities occur in the opposite case. 



We can again restrict ourselves to the case |n| = 1, and 
therefore we have to find the eigenvalues A of the matrix 
E7i\, In the relevant case for which u\\ = U22 = u (see 
below) this can be done analytically. We obtain: 



2\l + u±\ l + 2u 



A = 1± 



Thus, we have the following scenario: 
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Stable 
Energetically unstable 
Dynamically unstable. 



The eigenvalues as well as the stability diagram is shown 
in Fig. 2. For u > U12 and when the interaction energy 
becomes more important than the (rotational) kinetic en- 
ergy {u ^> 1), the solution is stable as long as u — u\2 ^ 1- 
Since both u and u\2 are proportional to the number of 
atoms N, we have that by increasing N one can com- 
pletely stabilize the spin monopole. An interesting situ- 
ation occurs for u — 1/12; in that case the matrix ETii 
is not simple, which implies that the perturbation only 
grows linearly with time. On the other hand, the dynam- 
ical instability occurring for u < u\2 also appears for the 
homogeneous stationary solution and therefore it simply 
corresponds to a phase separation of the two components 
(internal states), as could be expected. 
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FIG. 2. (a) Excitation frequencies as a function of u for 
M12 = 10; (b) Stability diagram 

In order to establish if the monopole continues being 
stable in a realistic 3-dimensional situation, we consider 
a trapping potential of the form [Fig. 1(c)] 



V(r, z) 



1 



2 2 
muj„z 



1 



muj 2 r 2 + Voe- r \ 



(9) 



This corresponds to a standard dipole trap with an off 
resonant Gaussian laser beam propagating along the z 
direction ||. Here Vq gives the ac-Stark shift at the 
center of the trap. The equilibrium point R along the 
radial direction is given by R — 2a 2 \n[Vo/(mu) 2 a 2 )]. 
In order to have a ring trap, we choose the parame- 
ters such that R ;§> Ar, where Ar is the typical size 
of the atomic cloud along the z and r directions. Writing 
\9(r,z,0)) = Mr,z,6)\ T) + ip 2 (r,z,0)\ |), the Gross- 
Pitaevskii Equation governing the evolution is 
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Here Uij = Anh 2 aij/m with a. 



being the s-wave 



scattering lengths corresponding to the different colli- 
sions, and 5i,2 two constant offsets. The functions ^>i, 2 
are normalized to the number of atoms N12 i n each in- 
ternal state. 

Let us first make the connection between the 3-D 
model and the ring. We consider the simple situation 
in which the motions along the radial and z direction are 
frozen. The conditions of validity of this situation will be 
discussed below. If the number of atoms in each internal 
level is the same N\ = N 2 = N/2, we can reduce the full 
three dimensional problem to the ring case studied above 
by writing ^ l<2 {r, z,9,t) = {N/ (Air)} 1 / 2 f h2 (r, z)0i, 2 (0, t), 
multiplying the coupled Gross-Pitaevskii equations by /* 
and integrating in r and z. We obtain that <px t 2 satisfy 
the Eqs. (0) with 



dz / rdr|/,(r,z)| 2 \f 3 (r,z)\ 2 , (10a) 
-00 Jo 



Oi = -0 (ndi +ei), t 
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t. 



(10b) 



where e% is the expectation value of the kinetic plus po- 
tential energy with the wavefunction /j. 

In general, we look for stationary solutions of the form 



Mr,z,9,t) = VN/2f 1 (r,z)e-^ t , 
Mr,z,9,t) = ^N/2f 2 (r,z)e lB e^ t , 



(11a) 
(lib) 



with /i,2 normalized real functions satisfying £1,2/1,2 = 
where 
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+ V(r, z) 
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We will concentrate in the case fin = U22 = 2 (equiv- 
alently, an = 022 = a) which is the relevant experi- 
mental situation described below. In the limit R 3> Ar 
the centrifugal term in L2 can be approximated by a 
constant Ti 2 /(2mR 2 ). If there is no phase separation 
(a > a 12 ) we find /i(r, z) =_/ 2 (r,j») = f{r,z) and 
fl = U2 — Ui = n/(2mR 2 ) + 82 — 5\. Solution ( pT| ) is again 
an eigenstate of ft -a with ft = [cos(0 — fit), sin(# — fit), 0], 
being the eigenvalue proportional to |/(r, z)\ 2 . Thus it 
represents a state where the atoms have their spin in the 
x-y plane forming a 2D monopole. The length of the 
spin at each point depends on the corresponding local 
density. We analyze now the stability of such a solu- 



tion. As before, we linearize around the solution (11) by 



adding a small quantity ai.2 and expanding it in powers 
of e ie . We obtain (§) where now H„ = L + K n + H lnt 
a 4 x 4 matrix operator with L = diag(Li, L±, L2, L 2 ), 
K n = h 2 /(2mr 2 )K n , and H int = f 2 H int .^ Using again 
R ^> Ar we have L\ — L2 and therefore L = L\ times 
the 4x4 identity matrix. Since we are only interested in 
the stability, we analyze the positivity of Ti. n . As before, 
we just have to study the cases \n\ < 2. Let us distinguish 
two situations: 

Weak interactions: in the limit N(a + &12)/ R <C 1 we 
have that the interaction energy Nu/[2irR(Ar) 2 ) is much 
smaller than the harmonic oscillator quantum fiuj. In 
that case, Ar ~ ao where ao = ^/(mu;)] 1 / 2 is the size of 
harmonic potential ground state. The spectrum of Ho,i 
is dominated by L. The radial and z dependence give 
rise to excitation energies khto (k integer). The lowest 
excitations k — correspond to a(r, z) oc f(r,z), and 
therefore we obtain A = h 2 /(2mR 2 )X (<C tko in absolute 
value ) , where A is given in (|8|) with u and M12 given in 
(|l0a| ) . Thus, for excitation energies lower than Hoj the 
problem fully reduces to the ring case. 

Strong interactions: In the opposite limit, we are in 
the Thomas-Fermi regime, where Ar = ao[32N(a + 
ai 2 )/-R] 1//4 - Now, one cannot simply separate radial and 
z excitations from ring excitations. The excitation spec- 
trum of Hq,i is dominated by L + H lnt . It is conve- 
nient to diagonalize H mt , and consider the cigcnfunc- 
tions separately, (a) Consider a^ n > — (51,32) <8> (1,-1): 
in this case H lnt is zero, and therefore the excitation fre- 
quencies correspond to those of L, which are of the or- 
der of kh 2 /[2m(Ar) 2 } > Ti 2 /{2mR 2 ). (b) For = 
(ff; 9) ® (1; l)i H mt gives (u + 1112)/ 2 , whereas for cjW = 
(g, —g) ® (1, 1) it gives (u — M12)/ 2 : the lowest energy will 
be of the order of N(u — ui2)/[2irR(Ar) 2 ]. As long as this 
energy is larger than h 2 / (2mR 2 ) we can consider sepa- 
rately the cases (a) and (b) treating K n as a perturba- 
tion; in both the correction is positive, i.e. the monopole 
is stable. In the opposite case, one has to be more care- 
ful in the perturbation analysis, since one cannot sepa- 
rate the cases (a) and (b); the excitation energies may 
become negative. Thus, we obtain a necessary condi- 
tion for stability N(u - ui 2 )/[27ri?(Ar) 2 ] > Ti 2 /(2mR 2 ). 
Using ( 10a ) we can write this condition as u — 1*12 > 1, 
which coincides with the basic stability condition derived 
for the ring. It means that the interactions have to be 
sufficiently strong to stabilize the monopole. 

In order to be specific, we will propose now a particular 
configuration to create the spin monopole. We consider 
an alkali atom in a ground F = 1 hyperfine state. We 
will assume that the energy of the mp = level is made 
higher (by using an off-resonant laser or radio-frequency 
field), so that it is not involved in the dynamics. In 
this case we can identify | f) = \F = l,mp = 1) and 
I I) = \F = l.mp = —1), the collisions do not change 
spin, and u\\ = U22 = u. Although the current status 



3 



does not allow to make a clear statement about the sign 
of u — u\2 it seems that for 23 Na it is positive, which 
is also agrees with recent experimental results concern- 
ing miscibility [H. In order to generate the 2-D spin 
monopole ( |TT| ) we propose to use an off-resonant Raman 
beam. The atoms are initially condensed in the internal 
| |) state. A Raman laser that connects the states | f) 
and | |) is then switched on. It should have the appro- 
priate spatial dependence so that the angular momentum 
in the z direction is changed by one unit jOJ. Denoting 
by fl(z,r,9) the effective Rabi frequency, the evolution 
equations are the above Gross-Pitaevskii Equations but 
with a coupling term between ipi and tp2 proportional 
to f2; the Raman detuning is incorporated to the defi- 
nitions of 6i t 2- Initially, one takes Si 5\, so that the 
laser does not affect the internal atomic state, since it is 
effectively out of resonance. Then, 62 ~8± is changed adi- 
abatically until 62 — 5\ — 0. The method is more robust 
than the one used to generate vortices since the spatial 
wavefunctions f 12(^1 z ) remain practically constant with 
our setup. Actually, in this case one can simply use a 7r/2 
laser pulse, taking 8\ = 82- This allows to generate the 
monopole state in a much faster time scale, which will be 
of the order of several inverse trap frequencies. 
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FIG. 3. Preparation of the 2D spin monopole. Trap 
parameters: lu z = 2u> r , R — 10ao, o = 5ao, 
Vo = 200hu! r , where ao = (Ti/muj) 1 / 2 ; interactions: 
fin = U22 = U12/O.9 = 36007kiVaoi Laser parameters: 
Q(r) = S7o[sin(fcx) + isin(ky)] with Qo ~ 0.28o; r and 
k = tv/(6R). (a) Final density distribution k/>i| 2 (dashed 
line) and \ip2\ 2 (solid line) at y — z — as a function of 
x; (b) Evolution of the population of the | f) (dashed line) 
and | J.) (solid line); (c) Evolution of the spin density after 
the preparation for taking 82 — deltai = 2u) r ; the triangles 
point along the expectation value of (a) and are proportional 
to the local density. 

In order to evaluate our proposal, we have performed 
a 3-dimensional numerical simulation of the Gross- 
Pitaevskii equations in the presence of the laser for the 
creation of the spin monopole. We have used an opti- 
mized three-dimensional collocation Fourier method with 
typically 80x80x40 collocation points and integrating in 
time with a symmetrized split-step operator technique. 



The results are shown in Fig. 3. In Fig. 3(a) we have 
plotted the final density profile along the x axis at the 
end of the generation process, whereas in Fig. 3(b) we 
have plotted the evolution of the population of the in- 
ternal levels until the monopole is generated. After this 
process, we switch of the laser and apply an internal en- 
ergy shift so that the spin start precessing, as is shown 
in Fig. 3(c). For this figure we have taken realistic pa- 
rameters. For example, taking Na with u r = 100Hz, 
ao = 2/im, and a — 52clb we have that the number of 
atoms is of the order of 2 x 10 5 . In order to ensure the 
stability, we have evolved the formed state in imaginary 
time (renormalizing the state after each evolution step) 
for about 20 trap oscillation times, without noticing any 
instability. 

In summary, we have shown that when the interactions 
are sufficiently strong, a 2D spin monopole becomes sta- 
ble in a ring trap. We have performed the stability anal- 
ysis, both in a ring situation as well as in the full 3D 
case. We have shown a way to prepare such a state, and 
verified with a full 3D numerical simulation that one can 
prepare it with current experimental parameters. 
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